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convenience of the numerical method, we introduced variable splitting rates repre- 
senting the details of the dynamics as model parameters which were not considered 
in previous works. The dependence on model parameters was studied, and it was 
found that the sea quark flavor asymmetry weakly depends on model parameters. 
It reflects the statistics principle contributes the dominant part of the asymmetry 
and the effect caused by details of the dynamics is small. We also applied the Monte 
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metries in kaons, octet baryons X, S, and A baryons, even in exotic pentaquark 
states. 
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SEA-QUARK FLAVOR ASYMMETRY FROM STATISTICAL BALANCE MODEL 

Although the proton is the simplest system in which the three colors of QCD neutralize 
into a colorless bound state, we still do not know how to describe the proton in terms of its 
fundamental quark and gluon degrees of freedom from basic principles. The structure of the 
proton is rather complicated due to the nonperturbative and relativistic nature of the quark 
and gluon in the protons. The complication also comes from the presence of sea quarks in 
the proton. The sea flavor symmetry naively assumed in the Gottfried sum rule |l( , which is 
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a symmetry between the light flavor u and d sea quarks inside the proton, was disproved by 
experiments of both deep inelastic scattering and Drell-Yan processes {2 — V | - 

Many theoretical attempts have been made to describe the origin of the nucleon sea and 



its antiquark asymmetry f?-l2l|. It is assumed that the primary mechanism to generate the 
sea is gluon splitting into uu and dd pairs. Field and Feynman 22] suggested that the extra 
valence u quark in the proton could lead to a suppression of g — >■ uu relative to g — >■ dd via 



Pauli blocking. But a subsequent calculation 



23] found that the effects of Pauli blocking 



are very small, and this result has been confirmed by another calculation [24|. Thus, it is 
believed that there must be a nonperturbative origin. For example, the meson-cloud inside 
the nucleon can account for such asymmetry [?-15] and chiral quark models 16-Ji]]. Also 
the large- A^ c approach can explain the flavor asymmetry of the antiquark distribution. 

Another attempt to understand the sea flavor asymmetry of the proton is from a pure 
statistical consideration in a kinetic equilibrium model 21| or "statistical balance model" as 
called in previous papers. The idea is rather simple and perspicuous: while the sea quark- 
antiquark uu and dd pairs can be produced by gluon splitting with equal probabilities, 
the time-reversal invariant processes of the annihilation of the antiquarks with their quark 
partners into gluons are not flavor symmetric due to the net excess of u quarks over d quarks. 
As a consequence, the u quarks have a larger probability to annihilate with the u quarks 
than that of the d quarks, and this brings an excess of d over u inside the proton. Taking 
the proton as an ensemble of a complete set of quark-gluon Fock states, and assuming the 
probability of 'arriving in' one state from others equals to the probability of 'leaving' it, 
one can obtain the probabilities of finding every Fock state (state density) in the proton. 
Thus one can calculate the quark and gluon content of the nucleon from a pure statistical 
consideration. It is interesting that the model gives a sea flavor u and d asymmetry as 
[d — u] ~ 0.132, which agrees with the experimental data 

The following diagram can describe the 'state shifting' between states. 
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Assuming kinetic equilibrium, we have these kinetic equilibrium equations: 

n n 

H c ijPi = H c iiPh (!) 

where pi is \i > state density, Qj is the non-normalized state-shift probability(NSSP) of 
\i >— > \j >,n is the total state number. Also there is the normalization condition 

n 

E* = l- (2) 

i 

If we know c^, we can derive state densities pi's by solving a system of n linear algebraic 
equations when n is a finite number. If n is infinite, we can get pi by asymptotic approach 
in some case if pi converges as n — > oo. Actually, if we change c^- to c^/Co, where C is a 
arbitrary constant, the result would be the same. It means we only need the ratios of NSSPs 

C{j s. 

If only considering the particle numbers of quark, anti-quark and gluon, the proton state 
can be described as an ensemble of Fock states 

\uud >, \uudg >, \uuduu >, \uuddd >, \uuddd >, • • • 
• • • I AT Nj N- Nj N > ■ ■ ■ 

Because the u-quark number N u = N u + 2, and N d = 1, all Fock state can be denoted 
with just three numbers as \N a , Nj, N g >. 

In order to derive the state density p\N a ,Nj,N g > we should know the probability of states 
shifting. We introduce the rate f q -+ qg as a quark splitting ability factor, there are 2N U + 
2Nj + 3 quarks(including antiquarks) in the initial state, so the NSSP of \N a , iVj, N g >— > 
\Nv,Nj,N g + l> is 

(2N u + 2N d +3)f q ^ qg . (3) 

We also introduce the splitting rate f g ^ qq and f g ^ gg , so the NSSP of \N n , Nj, N g >— > 
\N n , + l,N g -l> and \N U , Nj, N g >^ \N U + 1, N d , N g - 1 > is 

Ngfr+M ( 4 ) 
and the NSSP of \N a , N d , N g >-> \N a , N d , N g + 1 > is 

Ngfg^gg- (5) 
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Now, we consider the time-reversal process and assume those fusion rates 

fqg^q = fq^qgi 
fqq^g = fg^qq^ 
fgg^g = fg^gg 

for time-reversal invariance. 

Hence, the NSSP of \N C , N a , N g >^ \N a , N s , N g - 1 > is 

(2N U + 2N 3 + m 9 fqg^q + N9{N f l) fg 9 ^ (6) 
the NSSP of \Nu, N s , N g >^ \N U - 1, N 3 , N g + 1 > is 

(N Q + 2)N a f q ^ g , (7) 
the NSSP of \N a , N s , N g >^ \N U , N 3 - 1, N g + 1 > is 

{N 3 +l)N 3 f qg -^ g . (8) 

We can see that the probability of uu annihilation is larger than dd annihilation in all of 
the proton states because of valence quark asymmetry. This is the origin of the sea quark 
flavor asymmetry. 

It is assumed that all the splitting and fusion rates are the same in the previous papers 21] . 



If we get all the non-normalized state-shift probabilities c^-, the state densities can be derived 
out if the particle numbers N u3g are finite. We set an artificial limit N a3g < N max and 
solve the finite linear equations. The numeric state densities are then derived. The sea 
quark flavor asymmetry can be written as: 

[d-u] = J2( N d- N u)P\Nu,NaA g >- (9) 
u,d,g 

The sea quark flavor asymmetry converges to 0.133 when N m g X increases. The result is 
consistent with experiment data {2-6]. Some subsequent works [251, 13 followed the kinetic 
equilibrium principle to study the spin of nucleons and the parton distributions in the proton 
and pion, and obtained quite good results agreeing with the corresponding experimental 
values. 

However, in the previous works, we assumed that all the splitting rates are the same 
as fq^qg = fg^qq = fg->gg = 1 an d did not estimate the "error bound" caused by the 
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assumption. As we can imagine, if the splitting-rates vary in different orders of magnitude, 
the convergence of flavor asymmetry will be bad. It is necessary to solve large N max linear 
equations. So we need a convenient numerical method to explore the effects of different 
splitting-rates and to study more complex hadronic states. 

MONTE CARLO SIMULATION APPROACH OF A KINETIC EQUILIBRIUM 

ENSEMBLE 

Monte Carlo simulation also can give the numeric state densities instead of solving alge- 
braic equations even when the number of states is infinite. Here, we want to explain some 
details about the Monte Carlo evolution on kinetic equilibrium and prove the equivalence 
between Monte Carlo evolution approach and solving algebraic equations. Let us start with 
an arbitrary initial state \i >, and then let it make a possible shifting during each unit step. 
The probability of the state \i > shifting to \j > is Cij/Co. Here, Co is an arbitrary large 
constant we introduced to ensure that the total shifting probability for each prior state is 
less than 1. It is required that C > J2j^i c ij f° r an prior states \i >, so the probability of 
staying in the prior state \i > is 

l-E^/Co- (10) 

The state evolves step-by-step as random walk, and we record the number of iteration steps 
as Ti while the state \i > is emerging. And after a large number of iteration steps T, the 
normalized \i > emerging probability is Tj/T. For each step while the state is \i >, the next 
step has the probability c^/Co to be \j >. So there are the times TiCij/Co of state shifting 
\i >— > \j >. Of course, other states also can shift to \j >, meanwhile \j > has chance to 
stay at \j >. That means the number of those steps \j > emerging should be 

T 3 = £ cy/CM + (1 - E Cfi/Coft. (11) 

The equation can be reduced to 

E c ^ = E c ^r (12) 

The equation is independent of the constant C . The value of C only determines the number 
of iteration steps needed to arrive at the equilibrium state after starting from an arbitrary 
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initial state. We can find the above equation is just the kinetic equilibrium equation (pQ), if 
we consider that the normalized \i > emerging probability Tj/T is equivalent to the state 
density as 

T./T = ^ (13) 

And we also have the sum condition 

E T * = T > (14) 

i 

which is equal to the normalization condition Eq(j2J). Hence, we proved the equivalence of 
the Monte Carlo simulation approach and solving algebraic equations. 

The Monte Carlo simulation approach provides a powerful method for solving kinetic 
equilibrium ensemble problems. This method is error-controllable and very useful especially 
on complex multistate systems, such as the applications to other hadrons in the following 
sections. We gain the same value of the sea quark flavor asymmetry 0.132±0.02 in the proton 
as expected. Here, the error bar ±0.02 is the standard deviation of results with different 
random number series, and the deviation will decrease when computing time increases. 



DYNAMICS-NONSENSITIVE SEA QUARK FLAVOR ASYMMETRY IN 

PROTON 

The fusion rate should be the same as the splitting rates for a time-reversal process. In 
other words, the evolution in the proton should be time-reversal invariant. But there is no 
principle requires that the quark and gluon splitting evolution abilities of g — > qq(gg) and 
q qg are equal. Therefore we should introduce three splitting-rates f q ^ qg ,f g ^ qq and f g ^ gg 
to represent the quark and gluon splitting evolution abilities which are determined by the 
dynamics of quarks and gluons. Each rate enhances the corresponding splitting or fusion 
evolution probability. In previous works, we assumed that all the splitting-rates are the 
same to be f q -+ qg = fg^qq = fg^gg = 1 an d did not estimate the "error band" caused by the 
assumption. In the present work, we introduced a numerical Monte Carlo approach. This 
new method is easy to apply to complex systems, and it is easy to put the variable splitting 
rates in evolutions and calculate the deviation caused by them. 

In the above section, we can see that the state densities or results are independent of the 
constant C . The numerical value of f g ^ q q, for example, is input as f g ^ qq /C . Therefore 



7 



the result does not depend on the absolute value of f g ^ qq - It means that the sea quark 
asymmetry does not depend on the absolute values of those splitting rates. Only two ratios 
between three splitting rates will affect the state densities and the value of sea quark flavor 
asymmetry. So, we can fix the rate f q ^ qg = 1, and vary the other two ratios f g ^ qq / fq-> qg 
and fg-> gg / fq^qg as two parameters in the model. 



TABLE I: The values of sea quark asymmetry for different ratios of splitting rates 

[d-U] X 100 fg-+qql fq^qg 



fg^gg/ fq^qg 


100 


10 


1 


0.1 


0.01 


0.001 





123 ±2 


124 ±2 


124 ±2 


124 ±3 


125 ±3 


126 ±6 


1 


131 ± 2 


132 ±2 


132 ±2 


134 ±3 


135 ±3 


136 ±6 


2 


137 ± 2 


138 ±3 


140 ±3 


140 ±4 


141 ±3 


141 ±6 


5 


150 ±2 


152 ±3 


153 ±3 


154 ±3 


156 ±4 


156 ±7 


10 


161 ±3 


163 ±3 


164 ±4 


164 ±3 


165 ±5 


166 ±8 


100 


179 ±4 


180 ±4 


180 ±4 


180 ±3 


181 ± 5 


182 ±9 



In Table [U the values of sea quark asymmetry for different ratios of splitting rates are 
listed. The previous result 0.132 ±0.02 is reproduced when f g -+ q q/ f q -+qg = fg-*gg/ fg-^qg — 1- 

From Table [U we can see that the asymmetry value [d — u] is not sensitive to the model 
parameter f g ^ qq / f q ^ qg , it is almost fixed when f g ^ qq / f q ^ qg varies in a very large range 
over five order of magnitudes. We also can find that the values of asymmetry are always 
larger than 0.123 whatever the splitting rates vary in an arbitrary large range. It reflects 
the principle of statistics contributes the dominant part of sea quark flavor asymmetry. The 
asymmetry only has a variation [c? — tZ] = (0.12 — 0.16) which is within 30% when f g ^ gg / f q ->qg 
varies in the range < f g -> gg / f q -+qg — 10, and still a small variation [d — u] = (0.12 — 0.18) 
even when f g ^. gg /f q ^ qg varies in a larger magnitude range < f g ^ gg / f q ^ qg < 100. So the 
effect brought from details of the dynamics is small and within the bound of the experiments' 
uncertainty. 

By now, we do not consider the probability of g — > ggg splitting and ggg — > g recom- 
bination yet, because the probability is suppressed by coupling constant and "three-body" 
splitting kinematics, g — > ggg can be regarded as two successive g —> gg, and its effect 
is same as the effect of increasing f g -+ gg , as we can see from Table [TTJ However, the rate 
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of three-body splitting g — > ggg must be much smaller than two-body splitting g — > gg or 
q —> qg, because the three-body phase space in perturbative QCD is suppressed by a factor 
of 2-3 order of magnitudes comparing with the two-body splitting. Though the parton split- 
ting in hadrons is a strong-coupling non-perturbative process, we believe that we still can 
safely assume f g ^ ggg /f q ^ qg <C 0.1 which only causes a very small enhancement as shown in 
Table [TTJ The effect of the splitting g — > ggg is thus negligible. 

TABLE II: The values of sea quark asymmetry [d — u] x 100 for different value of f g ^ ggg / f q ^ qg , 
for fg^ qg = 1, fq-^gg = lj fg^qq = 1 an d fg^ggg = fggg^tg 



fg^gggl fg^ig 





0.1 


0.2 


0.4 


0.6 


0.8 


1.0 


132 ±2 


135 ± 2 


137 ± 2 


142 ±3 


145 ±3 


148 ±3 


150 ±4 



Because the effect of the splitting g — » ggg and recombination ggg — > g is negligible and 
the asymmetry value of [d — u] is almost independent of the parameter f g ^ qq /f q ^ qgi there is 
only one parameter fg^gg/fq^qg can vary the asymmetry. This parameter is QCD relevant 
and it is the only input from dynamics. If the parameter could be fixed by analysis of QCD, 
the deviation on sea quark asymmetry caused by the details of dynamics can be determined 
and the sea quark flavor asymmetry in proton is predictable. 

These two splitting vertices are QCD vertices and have the same coupling constant. The 
splitting kinematics of g — y gg and q — )■ qg are also similar. So, the splitting rates of g — > gg 
and q — >• qg should be in the same order of magnitude. The assumption can be supported 
by the integrations of Altarelli-Parisi(A-P) splitting functions. Though these equations are 
valid in the perturbative region and the parton splitting in hadrons is a nonperturbative 
process, the ratio of the total splitting rates is still inspirational. The ratio parameter 
fg-+ gg / fq^qg can be heuristically "derived" from Altarelli-Parisi splitting functions [26]. 

The A-P splitting functions are 

P(q q(z)g) = C F \±^, 
1 — z 

P{g g(z)g) = C A [^ + -^- + z (i- z )], 
P(g^q(z)q)=T R [z 2 + (l-zn 



where the color factors Cp = 4/3, Ca = 3 and T R = 1/2. 
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The integrations of A-P splitting functions are assumed to be the total probabilities of 
quarks and gluons splitting. So the splitting-rates directly to be: 



The rates f q -+ qg and f g -> gg are logarithmic divergent when the integration limit z min — > 0, 
but fortunately the ratio between the two rates is not divergent, and thus we have the model 
parameter 



when z min — > 0. The ratio parameter is not sensitive to the integration limit z m i n . For 
example, when z m i n = 0.1, the ratio is 2.01 which is close to 9/4. Such small deviation 
change on parameter f g ^ gg / f q -> qg dose not have effect on sea quark asymmetry. Considered 
the integration limit is relative to Q 2 scale, then the model parameter f g ^ gg / f q ^ qg and sea- 
quark asymmetry are not sensitive to Q 2 scale. We estimated the ratio parameter by the 
perturbative A-P splitting functions, it is just the ratio of color factors. We assume the 
parameter value is still similar in the nonperturbative region. 

The nonsensitive parameter f g ^ qq / f q -> qg also can be derived by above method. But, it is 
relevant to the integration limit or Q 2 scale. The dependence can be extracted as ~°j 075Tfl 



when z m i n is small on the order of magnitude and becomes zero when z m i n — > 0. For example, 
the value of parameter f g ^ qq /f q ^ qg = 0.005 when z min = 10~ 6 , and the value is not sensitive 
to the magnitude of z min or Q 2 scale because of its logz min dependence. We can see from 
TableJH the sea-quark asymmetry is not sensitive to this parameter even it is so small. 

As discussed above, the ratio f g ^ gg /f q ^ qg is almost fixed to ratio of color factors as 9/4 
and the asymmetry is independent of other details except the parameter f g ^ gg / f q ^ qg . There- 
fore we arrived at the following conclusion: after considering the detail of QCD especially 
the color factors, we can predict the sea quark flavor asymmetry in proton is 0.142 ± 0.03. 
It is enhanced a little compared to the value given in the previous papers. More precise 
measurement of [d — u] is needed to examine the statistical balance model. 





l-z., 



P{g -j. g(z)g)dz Ca = 9 
P{q q{z)g)dz C F 4' 



l-z. 
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The x-dependent d(x) — u(x) can be derived from deep inelastic scattering and Drell- 



Yan processes, and Jq d(x 



u[x\ 



to x — > and 



dx is given by extrapolating d(x) — u(x) 
x — > 1. The sea quark asymmetry values from three collaborations are listed in Table. IIII1 
they are all consistent with the sea quark asymmetry value predicted above. The value of 
E866 seems a little bit smaller compared to the prediction value, but the x range of the 
E866 measurement is narrow and the uncertainty brought by extrapolating to small x is out 
of control. So, more precise measurements are needed to test the prediction. 



TABLE III: J [d(x) — u(x)] dx as determined by three experiments. The range of the measurement 
is shown along with the value of the integral over all x (Q 2 = 54 GeV 2 /c 2 ). 



Experiment 


x ran 


ge 


Jo [ d ( x ) ~ u ( x )] dx 


E866 


0.015 < x 


< 0.35 


0.118 ±0.012 


NMC 


0.004 < x 


< 0.80 


0.148 ±0.039 


HERMES 


0.020 < x 


< 0.30 


0.16 ±0.03 



SEA QUARK FLAVOR ASYMMETRY IN MESONS 

Because the sea quark asymmetry value is not sensitive to details of dynamics and only 
depends on the parameter f g ^ gg / f q ^ qg which is almost fixed as 9/4, then it should not only 
work for the proton, but also for the mesons and other baryons. We suppose the statistical 
model also has validity on predicting sea quark asymmetry in other hadrons. M. Alberg , E. 



M. Henley 



and C.-B. Yang 



28] derived the parton distributions of pions according the 
statistical model, but the sea quark asymmetry is zero because of the same valence quark 
number in pions. While the valence quark numbers of the u and d quarks are different 
for the kaons, for example, K + (us) has one u valence quark and no d valence quark. The 
statistical balance model predicts the sea quark asymmetry value d — u = 0.284 in K + , 
when f g -> gg / f q ->qg = 9/4. In the same way, the sea quark asymmetry value [d — u) = —0.275 
in K°(ds) and [d — u] = —0.275 in K°(ds), d — u = 0.275 in K~(us). These sea quark 
asymmetry values are also not sensitive to dynamics as shown in Table IIVI 

We can see from Table HVl that the asymmetry [d — v] is independent of f g ^ q q/f q ^ qg and 
varies in a small range 0.263-0.31 as f g ^ gg / f q ^ qg varies in a large range 0-10. 
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TABLE IV: The values of sea quark asymmetry d — u in K + for different split factors 



[d - u] 


fg-*qq/ fq^qg 


f If 
Jg^ggl Jq^qg 


1 00 


1 


i 

_L 


1 


01 


n 


0.263 


0.264 


0.264 


0.264 


0.265 


0.1 


0.264 


0.265 


0.265 


0.266 


0.266 


1 


0.272 


0.274 


0.275 


0.277 


0.278 


5 


0.296 


0.300 


0.303 


0.304 


0.305 


10 


0.311 


0.312 


0.312 


0.312 


0.313 



SEA QUARK FLAVOR ASYMMETRY IN BARYONS 



We also use our statistical model to predict sea quark asymmetry for baryons. In a 
previous paper [29[, L. Shao et al. derived the octet baryons' sea quark asymmetry values 
by the method of solving linear equations. They give [d — u] = 0.41 in H + (uus) and 
[d — u] = 0.276 in S + (wss). In this paper, we get the same number by the Monte Carlo 
approach. We can find that the sea quark asymmetry value in S + (iiss) is almost the same 
as the meson K + (us) because their u and d valence quark numbers are the same. So, in 
the statistical model, the s valence quark number in the hadron has a negligible effect on 
the [d — u] sea quark asymmetry. We also find the sea quark asymmetry values in the octet 
baryons are not sensitive to details of dynamics, they just depend on the valence quark 
numbers in those baryons. The asymmetries [d — u] in Y, + (uus) and 3 + (uss) are enhanced 
a little to be 0.42 and 0.285 when f g _> gg /f q ^ qg = 9/4. 

Besides octet baryons, we also derived A baryons' sea quark asymmetry value as: 



d-u = 0.50 for A ++ (uuu) , 
d-u = 0.14 for A+(uud) , 
d-u= -0.14 for A°(udd) , 
d-u = -0.50 for A~(ddd) , 



where, f g ^ gg / f q ^ qg = 9/4. The sea quark asymmetry in A + (uud) is the same as in 
proton because of their same u and d valence quark numbers. Of course, the asymmetry in 
A°(udd) is the same as in neutron. 
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We also derived exotic baryons' (pentaquark states) sea quark asymmetry values as : 

d-u = -0.14 for $ (ssddu), 
d — u — 0.14 for $~(ssuud), 

where, the sea quark asymmetry values are the same as in the proton because of their same 
u(u) and d(d) valence quark numbers. 

If there is such a pentaquark state X ++ (uuuds), then its sea quark asymmetry value 
would be [d — u] = 0.21 derived by the statistical model. 



CONCLUSIONS 



In the previous works in the statistical balance model, the sea quark flavor asymmetry 
[d — u] = / dx(d(x)—u(x)) in the proton was computed using the "linear equations method" . 
Because of the difficulty and limit of the linear equations method, it is hard to apply the 
method to more complex systems. It is also assumed that all the splitting-rates are the 
same, j q -+ qa = f g ^ qq = fg^gg = 1 in the previous works, and the "error band" caused by the 
assumption was not estimated. In the present work, we introduced a numerical Monte Carlo 
approach. This new method is easy to apply to complex systems, such as other mesons and 
baryons. We also introduced the variable splitting rates representing details of the dynamics 
and studied the dependence on them. We find the sea quark flavor asymmetry in the proton 
is always larger than 0.123 whatever the splitting rates vary over an arbitrary large range. 
It reflects the statistics principle contributes the dominant part of the asymmetry. The 
asymmetry is almost independent of the model parameter f g ^ qq / f q ^ qg and only changes 
within 30% when f g ^ gg / f q -+ qg varies in the range — 10. So the effect caused by details of 
the dynamics is small and within the bound of the experiments' uncertainty. However, these 
two splitting vertices are QCD vertices and have the same coupling constant. The splitting 
kinematics of g — >■ gg and q — » qg are also similar. So the splitting rates of g — > gg and 
q — >■ qg should be in the same order of magnitude. The assumption can be supported by 
the integrations of Altarelli-Parisi splitting functions. Though these equations are valid in 
the perturbative region, one may heuristically assume that the ratio of the total splitting 
rates obtained from them holds approximately also in the nonperturbative regime. The 
parameter f g -^ gg / f q ^ qg can be fixed to the ratio of color factors as 9/4 by integrations of 
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Altarelli-Parisi splitting functions. According to the above reasons, we can conclude that 
the prediction only from a statistics principle has an accuracy < 30%. Or, in other words, 
the details of the dynamics only bring less than 30% effect. After considering the details of 
QCD especially the color factors, the sea quark flavor asymmetry in proton is enhanced to 
0.142 ±0.03 which is consistent with present experimental measurements and can be tested 
by more precise measurements. 

The sea quark asymmetries are not sensitively dependent on the details of dynamics. The 
sea-quark flavor asymmetry derived only from statistic principle contributes the dominant 
part of the asymmetry. It strongly implies that the origin of the sea-quark flavor asymmetry 
of hadrons is the asymmetry of valence quarks. We also applied this Monte Carlo approach 
of statistical model to predict the sea quark asymmetries in kaons, octet baryons S, S, and 
A baryons, even in exotic pentaquark states. All these asymmetries just only depend on 
the valence quarks number in those hadrons. The sea-quark asymmetries for different u and 
d valence quark numbers are listed in Table |VJ These values can confirm the mechanism 
we proposed to explain the sea quark asymmetry in proton. It can be observed from Table 
El that the sea quark asymmetries are enhanced by the difference of corresponding valence 
quark numbers and suppressed by the sum of valence quark numbers. When the valence 
quark numbers [u v ] > [d v ] , the sea-quarks u are easier to annihilated because of the existence 
of more u valence quarks and it leads the sea quark asymmetry. On the other hand, the 
larger total number of valence quark [u v + d v ] suppresses the relative difference of valence 
quarks and weakens the sea quark asymmetries even if [u v — d v ] remains the same. These sea 
quark asymmetries for hadrons, except the proton, are listed purely for theoretical interest, 
as it is not known presently how to access this information in experiment. 
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